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Abstract 

We consider a brownian motion on a general graph, that starts at time t = from some 
vertex O and stops at time t somewhere on the graph. Denoting by g the last time when O 
is reached, we establish a simple expression for the Laplace Transform, L, of the probability 
density of g. We discuss this result for some special graphs like star, ring, tree or square 
lattice. Finally, we show that L can also be expressed in terms of primitive orbits when, for 
any vertex, all the exit probabilities are equal. 



Around the years 1930, P Levy got several arc-sine laws concerning the ID Brownian Motion 
on an infinite line. Let us consider such a process starting at t = from the origin O and 
stopping at time t and denote by g the last time when O is reached. The second arc-sine law 
discovered by Levy concerns the probability law of g. It can be stated as follows: 

2 [u 
P(g < u) = — arcsinw— (1) 

with the probability density Vt(u) = — = , (2) 

an 7r y/u(t - u) 

In particular, the Laplace transform of Vt{u) is written: 

roo rt 1 

/ dte-W dnV t (n)e^ u = (3) 
Jo Jo V717 + 6 



Levy also proved that the same law occurs when g represents the time spent in the region (x > 0) 
(first arc-sine law). 

The infinite line with only one point, O, specified, can be viewed as a kind of very simple graph 
consisting in one vertex, O, and two semi-infinite lines originating from O. 

Our goal in this Letter is to get an analogue of the second arc-sine law © but for a quite 
general graph. 



Let us recall that graphs properties have interested for many years physicists 2 a as well as 
mathematicians [3]. For instance, the knowledge of spectral properties of the Laplacian operator 
is useful to understand superconducting networks t 4. as well as organic molecules jSj. We also 
know that spectral determinants play a central role in the computation of weak localization 
corrections (Hj. Recently, graphs have also been used in the context of nonequilibrium statistical 
physics (see j7] and references therein) and also in the field of quantum chaos [S] . 

The study of Brownian motion on graphs is also an active area of research. For instance, 
occupation times E3 or local time distributions have been recently computed for general 
graphs. 



So, let us begin by setting the notations that will be used throughout this Letter. We consider 
a graph made of V vertices, numbered from to V — 1, linked by B bonds of finite lengths. The 
coordination number of vertex a is m a (X^=o 777-0 = 2-B). 

On each bond [a/3], of length l a p, we define the coordinate x a p that runs from (vertex a) to l a p 
(vertex /?). Sometimes we will also use the coordinate xp a = l a p — x a p. The set of coordinates 
{x a @} is simply noted x. 

An arc (ct/3) is defined as the oriented bond from a to ft Each bond [a(5] is therefore associated 
with two arcs (af3) and ((3a). In the sequel, we will consider the following ordering of the 2B 
arcs: (0/zi)(0/i 2 ) • • • (0/i mo ) . . . (aft) . . . (aft) . . . (a/3 m J . . . 

Now, let us consider a brownian motion starting at some vertex O (label 0) of the graph at time 
t = and stopped at time t. Let g be the last time when O is visited. The probability, Pt(g < u) 
can be written: 

Pt(g < u) = / dx' dx Gq(0, 0; x, u) Gd(x, u\ x' , t) (4) 

J Graph 

where Go is the free propagator on the graph and Ge> is the propagator with Dirichlet boundary 
conditions in O. Q simply expresses the fact that the brownian particle will not hit O between 

the times u and t. Let us define the probability density Vt{u) = — — ^t— 1 — -. In the sequel, we 

du 

will be especially interested in the Laplace transforms: 

/■oo r-t 

C = / die" 7 * / du P t (g < u) e~Z u (5) 
Jo Jo 

poo rt 1 

L = / dte-^ duVt(u)e^ u = + (6) 

Jo Jo 7 + ? 

With (jlj), we see that: 

C = [ dx'dx(o 1 x)(x — — — x'\ = [ dx'dx G (0,x;j + G D (x,x';j) 

JGraph \ ii+7 + ^ / \ i/fl + 7 / JGraph 

(7) 



1 d 2 

On the bond aft, H and Hd are equal to , and we define: 

2 dx z 



a/3 

r (af3) v r (n i C \ r M) r dg (0, x Q/3 ; 7 + 

Q « ^o(0,^;7 + ; Go =^m o ^— (8) 

g^\x)= } ^g D (x,x' a ^) ; g^(x)= hm n dgg( ;;^ ;7) (9) 

cup ap ap 
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The integration over x 1 in Q is easily performed with the help of the equation 
We get 

„ -1-1 mo 

/ dx' g D (x,x>; 7 ) = - - — £s D (0w) (*) (ii) 

•/Graph 7 2 1 ~[ 

To go further, we have to specify the behaviours of the resolvants Go and Qd in the neighbourhood 
of all the vertices. 

Let us consider some vertex a with its nearest neighbours i = 1,2, ... m a , on the graph. 
Suppose that the brownian particle reaches a. It will come out towards /3j with some exit 
probability p a ft- This implies the following equations to be satisfied by Qq: 



= fa , Va (12) 



y y _ y 



Pafr Paf3 2 Pa/3 ma 

Remark that the resolvant will be continuous in vertex a only when the particle exits from a 
with the same probability in all the directions. 

Moreover, current conservation implies: 

m a 



Y,Qo {aPi) = , if a^O (13) 

i=i 

*£g' (°rt = -2 (14) 



i=i 



On the link [a(5], (7o(0, x a p; 7 + <!;) must satisfy: 

/id 2 



2 rlr 2 

z ax aj3 



+ 7 + ? K?o(0,x a/3 ;7 + = (15) 



with the solution: 



° ' ™ ' ° sinh \/2(7 + 0^/3 ~" sinh ^2(7 + 

So, we deduce: 

Oo (aP) = ^2(7 + (-^(7 + + ^(7 + # Q) ) (17) 

with: 

Ca^(7 + = coth ^2(7 + 0^/3 = ^(7 + (18) 

^(7 + = — — =L=— = ^ a (7 + (19) 
sinh a/2 (7 + £) 'a/3 

Those considerations allow to write eas. (|121 ITlfl I14|) in a matrix form: 
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+ =^WW) 1 (20) 

where / and I are two (V x 1) vectors. The components of / are the quantities f a and for the 
components of I we have lj = 25 m. 

M("f + £) is a (V x V) vertex matrix with the non vanishing elements (a runs from to V — 1): 



M( 7 + 0*a = $>aftC aft ( 7 + £) ( 21 ) 
i=l 

Af(7 + 0«/3 = -P/9a«^(7 + if [a/3] is a bond (22) 

= otherwise (23) 

Similar considerations hold for Qd(x, x'; 7) and we get: 

# l} (x) _ _ s{f "^(x) 



Pa/3i Pa/3r, 



F a (x) , if a / (24) 



G ( D m) (x) = ; » = !,..., mo (25) 



gg^)(x) = , if (26) 
i=l 



and the matrix equation Mb{i)F{x) = — Id{x) (27) 

7 

Md is the matrix M with the first line and first column deleted. F{x) and Id{x) are (V — 1) x 1 
vectors. The components of F(x) are the quantities F a (x), a = 1, . . . , V — 1. The nonvanishing 
components of £d(x) are defined as follows. 

If x 6 [aft] with a/0 and ft ^ then 



(ir>(«o6))o = s ab(l) sinh(v/27 x fea ) (28) 

{lD{xab)) b = s a6 (7) sinh(y / 27 x a6 ) (29) 

If x € [0/ij] then 

(M 2 ^) = «o w (7) sinh(v^7 ^o w ) (30) 
Performing the integration over x in Q, we finally get the following simple expression for L: 

r_ 1 (^(7 + 0)oo 

v/TTT+O (M- 1 (7))oo 1 j 



This is the central result of this paper. 
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Remark that, for the special case when Qq and Gd are continuous Green's functions, we can 
recast (|31|) in the form 1 

L= lgo(0 0; 7 + {) 
7 6o(0,0; T ) 



To go further and deduce from (|31|) some general properties of the density Vt (u) , let us consider 
the following identity 

/ dt e~~ ft / du e^ u r{u)s(t - u) = TZ^y + £)«S(t) (33) 
Jo Jo 

where 1Z and S are the Laplace Transforms of r and s. 

(M~ 1 (n)) 1 
Setting K(p) = )£±>M an d s{p) = ^-^^-^ in (ESI), we realize from (EB and © 

Vp Vp ( M (P))oo 

that the probability density factorizes in the following way: 

Vt{u) = r{u).s(t-u) (34) 

This result holds for any graph. 

It is now easy to analyze the limiting behaviours of Vtiu) when u — > + and when u — > t~ . 

In the first case, we can write Vt(u) ~ r(u)s(t). Moreover, for large p values, M(p) — ► 1 and 

ft(p) „ _L. So 

Vt(u) = s(t) , when u -> 0+ (35) 



On the other hand, when u — > t , VAu) ~ r(t)s(t — u). For large p values, <S(p) ~ — and 
finally 

"Pt(u) ~ = r(t) , when u — > (36) 

V vr (t - w) 



Let us now develop some examples. 
For the graph of Figure 1, we get: 



vW E"=i po, th( 72(7+0^) 

when Zoi — ► oo , i = 1, . . . , n (38) 



lr rhis is readily done with the help of the Dirichlet Green's function: 

Qd[x,x 7) = y (a:,ai |7) 

C>o(0,0;7) 

Using completeness and also the relationship 



/ da;' Qo(x,x; X) = ~ 

J Graph A 



' Graph 

all the integrations in Q can be performed, leading to [1M2ti . 
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Figure 1: a star-graph with n legs of lengths loi, i = 1, . . . n, originating from O. 



When all the legs become infinite, L does not depend on the exit probabilities from O, poi, 
because in that case all the legs are equivalent. The same behaviour occurs for a finite symmetric 
star (Iqi = . . . = Iq„ = I) where we get 



th(v^Z) 



^7(7 + 6 th(^2(7 + O0 

that leads to the density 

+CO „ 2,2 \ / +OO 



(39) 



v t (u) = ; 1 y: (-i) p e-^ ] ( J2 ] (4°) 



j / (2m — 1) 7r (t — u) » ' 



— [ ^ e 5" ) ( J2 e"" 53 " J (41) 



\n=— 00 



This last formula allows for studying the large u and t behaviours of Vt(u). For instance: 

V t (u) ~ J2 when a » I 2 and (t - u) > I 2 (42) 

In Figure 2, we consider a ring where the vertices O and 1 are linked by two bonds of lengths 
1$ and 1$ (I = 1$ + Iqi). The exit probabilities from O (resp. 1) are p^ and p^ (resp. 



and p^q) and we define Cqi (7) = coth -^2jl$, Sm (7) = pr, « = 1,2. 

sinh x /2^/ V 




Figure 2: a ring with two vertices, and 1, linked by two bonds (1) and (2). 



Adding intermediate vertices on each bond, we easily show that the matrix M writes: 



6 



/ Z^i=i Pqi c i 



M 



01 



01 



v^2 («) («) 
Z^i=lPl0 S 01 



v^2 («) («) 
Z^i=lPl0 c 01 



When = = 1/2, L does not depend on pjg and we get: 



th(V27 



V / 7(7 + 6th(V2(^+0|) 



(43) 



In that case, the ring can be viewed as a star with two legs of equal lengths, 1/2, glued at their 
endpoints. Thus, (|4*3*|) is simply (|3T?|) where I has been replaced by 1/2. 



Remark, however, that we should obtain a quite different result when taking p^ = p 1 ^ = 1/2 : 

1 coth(V2( 7 + 0!) + ^(7 + £) 



,( 2 ) 



L 



Vt(7 + 6 coth(V2^4) + F{i) 



with 



F( 7 ) 



(rf )-pg)sinhV2^(e-«) 



cosh -v/27/ — 1 



(44) 



(45) 



Let us now turn to a brownian motion starting at the root O of a regular tree of coordination 
number z and depth n (see Figure 3 part a) for z = 3 and n = 2). We suppose that all the links 
have the same length / and also that 



Pa/3i = = — for any vertex inside the graph 



z 
1 



otherwise 



(46) 
(47) 




b) 



O 



Figure 3: part a) a tree with coordination number z = 3 and depth n = 2; part b) the equivalent 
linear graph with exit probabilities given in ()48|) and 1)49(1. 
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For the last passage 
probabilities 



problem in O 



this tree is equivalent to a (n + 1) linear graph with the exit 



POX = 1 = Vn.n-X (48) 

Pi,i-x = ~ , Pi,i+i = 1 - - , i = 1, • • • , n - 1 (49) 

z z 

With c( 7 ) = coth \f 7 2r{l and 5(7) = — - — 7=7, we can write the relations 

sinh y 27/ 

detM( 7 ) = c( 7 ) detM n ( 7 )-i(s( 7 )) 2 detM n _i( 7 ) (50) 

detM n _ fc ( 7 ) = c( 7 ) detM n _ & _ 1 ( 7 ) + - (- - lj (s( 7 )) 2 detM n _ fc _ 2 ( 7 ) , < k < n — (SI) 

detM 2 ( 7 ) = 1 + I( S ( 7 )) 2 (52) 
detMi( 7 ) = 0(7) (53) 

M n _fc is the matrix M where the (fc+l) first lines and first columns have been deleted. Recursion 

relations can also be written for (M _1 ) = -. In the limit n — > 00, we get: 

V J 00 detM 



c(7) (!-!) + ^/l+(l-|) 2 (K7)) 2 



L = -^=^= , (54) 

=(7 + (1 - 1) + ^1 + (1 - 1 f «7 + or 

As expected, z = 2 leads to the second arc-sine Levy's law. 

As a final example, let us consider an infinite 2D square lattice of stepsize / with exit probabilities 
equal to 1/4 in any vertex. A standard computation based on the tight binding model 12\ leads 
to: 

K 



K(\) is a complete elliptic integral of first kind [T3*] . 



To complete this work, let us come back to a general graph and choose equal exit probabilities 
in each vertex: 

1 , . 

Pa/3, = — ; 1 = 1, ...,m a (56) 



V a 



In that case, L can be expressed in terms of primitive orbits on the graph. Recall that an orbit 
C is said to be primitive if it cannot be decomposed as a repetition of any smaller orbit. From 
[Tl| . we get: 



n 

c 



1 



n 

[a/3] 



1 



detM (7) 



(57) 



The product on the left hand-side is taken over all the primitive orbits C of the graph. /(C) is 
the length of C and the weight //(C) is the product of all the reflection - or transmission - factors 
along C. (An orbit is a closed succession of arcs . . . (to) (a/3) . . . and a reflection (transmission) 
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occurs in a if r = f3 ( r 7^ (3). The corresponding factors are ( 1 ) for a reflection and 



2 \ 

) for a transmission). 

Now, imposing Dirichlet boundary conditions in O, we get, still with |14j : 



n(l-/i (0) (C)e-^ 5 )) 

c 



JJ (l - e-V^a/, 
[aj9] 



detA/bfr) (58) 



fj,(°\C) is the same as /x(C) except for orbits visiting O where the reflection and transmission 
factors are respectively — 1 and 0. 

It is now easy to realize that the ratio — ~t~t' only involves the primitive orbits Co that visit 

detM(7j 

O at least once. Thus: 



(l-M (0) (Co)e-^<°°> 
(l - KCo)e 



(M- 1 ^)) = - (59) 



Finally: 

r _ 1 rr / j ~ M (0) (^je-^+^^) \ n / 1 - MCo) e -y^(^ , (6Q) 
V7(7 + ~ \ l- M (0)(C )e-V^(O») / ~ \l-//(Cb)e-V^K)'(Co) 1 

j^J means that we only consider primitive orbits with reflections in O. 



As an example, let us see how (|6Uj) works with the star-graph of Figure 1 with n = 2 legs. 
Only one primitive orbit, with length 2(Zoi + Z02) and two reflection factors equal to 1, contributes 
to the product TT. 

Co 

Concerning two orbits, with lengths 21q\ and 21q2, will contribute. For each one, the two 
Co 

reflection factors are (—1) and 1. For L, we get the result: 

1 (1 _|_ e -V 2 (T+£) J M)(l + e -V 2 (7+e)«02N M _ e -2*S^i(l0l+l02)\ 

L ~ \/7v7 + T) (1 + e - 2 v^7'oi)(l + e -2V2^ 2) ^ _ e -av^+0(I i+laa)) ^ 
_ 1 cosh y2(7 + g)/ i cosh ^2(7 + ^)^02 sinh y^^pi + Z 02 ) 

^7(7 + cosh ^27/01 cosh ^27/02 sinh ^2(7 + Cj(?oi + I02) 

in agreement with 1)37(1 . 

For the ring of Figure 2, two (time-reversed) orbits of length I should contribute to Y[ an d only 

Co 

one of length 21 should contribute to TT. So: 

Co 



- 1 (64) 



v/7(7 + f) thV2( 7 + ?)l 
as given by (fHH) . 
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